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Abstract 

We construct the non uniform AdS black string solution with a perturbation theory in six dimensions, 
focusing on the backreacting second order correction. The backreactions at second order give the first 
relevant corrections to the thermodynamical quantities. Our results show that for configurations with 
horizon radius and length in the extradimension small compared to the AdS radius, the properties of the 
non uniform black string are similar to the locally asymptotically flat case. For black strings with small 
horizon radial coordinate but large length in the extradimension, the thermodynamical properties of the 
solutions are affected by the AdS curvature. 

1 Introduction 

The last decade has witnessed a growing interest for solutions in AdS space, largely motivated by the 
AdS/CFT correspondence [H[5]. This conjecture states that a solution of Einstein field equations which 
approaches asymptotically the AdS space is dual to a conformal field theory living on the conformal boundary 
of the AdS. In this context, black objects are of great interest; for example, the Hawking-Page phase transition 
PQ between the five dimensional spherically symmetric black holes and the thermal AdS background was 
interpreted by Witten, through AdS/CFT, as a thermal phase transition from a confining to a deconfining 
phase in the dual four dimensional JV = 4 super Yang-Mills theory [3] . 

In four dimensional asymptotically flat space, there is a uniqueness theorem on black holes, implying 
that the horizon topology of black objects is £2. In higher dimensions, this theorem does not hold and many 
efforts have been done in order to construct different kinds of black objects, such as black strings [TO] of 
horizon topology Sd-3 x Si, black rings [IT] (toroidal horizon topology), black holes (Sd-3) [9] and their 
generalisation with charge, rotation, cosmological constant, etc. 

Gregory and Laflamme [5] showed in '93 that black strings are unstable towards long wavelength pertur- 
bations in asymptotically flat space (see [7] for a review) . After that discovery, many efforts have been done 
to understand the end point of the black string instability in asymptotically locally flat spacctimcs. It was 
believed that the endpoint of the instability should be a caged black hole, where the asymptotical space in 
d + 1 dimensions is Md x Si (with Md the d-dimensional Minkowski space), but it was argued [5] that a 
transition between a black string and a caged black hole would take an infinite proper time at the horizon. 
This lead to the study of non uniform black string, a black object which is non translational invariant along 
the Si coordinate, first in a perturbative approach [14] (see also |15j ) and after in the full non linear regime 
|19j . The resulting thermodynamical phase diagram |18j is now well known although dynamical phase tran- 
sitions are still to be constructed. In this diagram, the non uniform and the caged black hole branch meet 
at the merger point, where a topological phase transition is expected to occur. 

More recently, it has been shown that this instability persists in asymptotically locally AdS space for 
small AdS black strings [13j where the ratio r^jl « 1, with the horizon radius and t the AdS radius, 
fixed by the cosmological constant A. There exists another phase of AdS black strings, namely large AdS 
black strings, where r^/i > 1 which is thermodynamically and dynamically stable [13] . The picture for these 
two phases is essentially similar to AdS black hole: small AdS black strings thermodynamically unstable 
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while large black strings are stable. As stated in [T], AdS space acts like a confining box: when the ratio 
rhjl becomes larger than some critical values, the wavelength of the instable mode cannot fit the 'AdS box' 
and thus all acceptable wavelength do not lead to instabilities. 

However, in asymptotically locally AdS spaces, much less is known about the counterparts of the solutions 
of A = Kaluza-Klein black objects The evolution of AdS black holes has been studied [TO] but still has to 
be done for AdS black strings. In order to construct a phase diagram in AdS, one has to consider various 
stationary solutions and compare their thermodynamical properties since these stationary solutions should 
be equilibrium configuration. This could give some clues on the endpoint of dynamical evolution and phase 
transitions between different black objects in AdS. 

In this paper, we consider the case of perturbative non uniform black strings in AdS6 (sometimes referred 
as pNUBS). We choosed to work in six dimensions because it turns out to be numerically easier than in 
odd dimension (log terms arise in odd dimensions and are more difficult to deal with) and in more than six 
dimensions (where the field decay is such that it is much difficult to have an accurate extraction of numerical 
quantities). The phase of non uniform black strings is connected to the uniform black string phase, since 
one can smoothly deform a uniform black string into a non uniform one. The perturbative approach will 
give the direction of the evolution of the thermodynamical quantities. 

This paper is organised as follow: In section [21 we present and explain the strategy for the perturbative 
approach. Section [3] is devoted to the asymptotic solutions for the background, the first order perturbation 
and the back reaction. We explain in details the numerical technique and the appropriate boundary condi- 
tions in section [4] and give the properties of the solutions in section [5] Finally, we expose and summarise 
our results in section [6] and [7l 



2 The model and the equations 

We consider the Einstein-Hilbert action with a negative cosmological constant in six dimensions, supple- 
mented with Hawking-Gibbons boundary term: 



where G is the six dimensional Newton constant, I the AdS radius, R the scalar curvature on the manifold 
M., g the determinant of the metric on M, h the determinant of the induced metric on the boundary dM. 
and K the extrinsic curvature of the boundary. 

In order to construct black string solutions, we consider the following ansatz for the metric: 

ds 2 = -b(r)e 2A ^dt 2 + e 2B ^ ( ^- + a(r)dz 2 ) + r 2 e 2C ^dn 2 . (2) 

V/M / 

We furthermore assume the extradimension z to be of finite size, z € [0, L] and the solution to be periodic 
in z. The ansatz of [12] corresponds to A = B = C — and leads to the uniform AdS black strings 
parameterized by the functions a, b, / which can be obtained numerically. 

With the ansatz ([2]), the Einstein equations reduce to a coupled system of partial differential equations. 
Since it might turn out difficult to solve the equations, we approach the full nonlinear problem with pertur- 
bation theory. To do so, we consider A ,B, C to be z-dependent perturbations of the uniform AdS black 
string. Accordingly, we parameterize the functions according to: 

A(r,z) = eA 1 (r)cos(kz) + e 2 [Ao{r)+A 2 {r)cos{2kz)} + 0(e 3 ) : 

B(r,z) = eB 1 {r)cos(kz) + e 2 [B (r) + B 2 (r)cos(2kz)} + O(e 3 ), (3) 
C{r,z) = eC 1 {r)cos(kz) + e 2 [C {r) + C 2 (r)cos(2kz)] + O(e 3 ), 

with e << 1 and k = 2it/L E W.. This small parameter e is related to the deformation parameter of the 
string, as we will see later. 
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The fields a, 6, f are the background fields and a solution for the background has been constructed in 
[T2] . The fields A\ , B\ , C\ are the linear perturbations and give access to the stability problem (T3] . The 
fields Aq,Bq, Co are back reactions, while A 2 ,B2, C2 are massive modes (see [Tl] for a construction of these 
fields in asymptotically flat space). The equations at second order in e decouple for the massive mode and 
the back reaction, the massive mode equations come with a cos factor while the back reaction equations do 
not. 

The equations for the background and first order can be found respectively in |12j and |13j . Let us remind 
the first order equations {B\ can be expressed in terms of A±, C\): 



A" = ol x A x + a 2 A[ + a 3 d + a 4 C[, C" = ip 1 A 1 + ip 2 A[ + tp 3 d + tp 4 C[, 



(4) 
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This is an eigenvalue problem, where k plays the role of the eigenvalue. In fact, the value obtained for k by 
solving the above equations is the critical Gregory-Laflamme wavenumber for wich the perturbation is static 
and fixes static perturbative non uniform black string length L since L = 2w/k. 

The equations for the higher order correction follows the pattern of those in asymptotically flat spacetime 
[M] . It is possible to solve Bo in terms of A), Co and of the lower order fields; the equations for the back 
reacting fields read: 
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where Ja , Jc are source terms. Their expression contain about seventy terms so we don't give them here. 
Let us just mention that Ja i Jc are quadratic in A±, C\ and depend on the background functions. 
Since equations ((5]) depend only on the derivatives of Aq and Co, the solutions are invariant under 



Aq — > Aq + const. , Co — > Co + const. 



(0) 



3 Near-horizon and asymptotics expansions 

The near-horizon behaviour for the background and for the first order is given in [12] (resp. |13j). we will 
remind them. We also remind the asymptotic expansion for the background fields and compute it for the 
first order. 

The near-horizon expansion for the background is given by [12] : 

a(r) = ao+a 10r/l ( r -r h )+0(r-r h ) 2 , b(r) = h (r-r h )+0{r-r h f , f(r) = ^ + ( r -r h )+0(r-r h ) 2 , 
(5r^ + £ 2 ) r h£ 

(7) 
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where ao, b\ arc numerical factors which are chosen in order to fit the following asymptotic expansion: 
o(r) = r ^ + l + l( l -) 2 + c A 3 +°( £ / r4 )> b (r) = ^ + l + l( e -) 2 +c t ( e -f + 0(£/r% 

/(r) = i + \ + l& + ^+ c ^ i+ °w r ^ ^ 

where c t , c 2 arc real constants to be determined numerically. 

The near-horizon expansion for the first order is of the form |13j : 

A 1 (r)=A 10 + A n (r-r h ) + O(r-r h f , Ci(r) = C 10 + C n (r - r h ) + 0(r - r h f (9) 



where An, C\\ are functions of A\§, C\q while the ratio A10/C10 must be determined numerically. Either 
A\o or Cio can be chosen arbitrarily since the equations are linear. 

Using the asymptotic form ([8]) of the background, the first order fields have the following asymptotic 
expansions: 
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where ip is a real constant determining the arbitrary normalisation of the fields. 

Before considering the asymptotic solution for the backreacting fields, we give their near-horizon be- 
haviour, since this expansion will be useful later: 

A (r) = Ao + A 01 (r-r h ) + ^(r-r h ) 2 +0(r-r h ) 3 , C (r) = C oa + C i(r^r h ) + ^(r-r h ) 2 +0(r~r h ) 3 , 

(11) 

where A$o, Coo, On are arbitrary real constants while 

2 (4 + 10^1) (A 2 10 + 6A 10 C W - 3C 2 10 ) k 2 r h 
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The functions A02, C02 are rather cumbersome and we refrain to give them explicitely. Here, we posed 
a Q = a(r h ), A w = Al(r h ), C w = Ci(r h ), C i = C'(r h ). 

The most general asymptotic expansion for the back reaction fields is given by: 
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with ao, ¥>o are arbitrary real constants and %, cr are real constants to be determined numerically 

Properties of the solution will be discussed in section [5] Some of these properties depends on the 
asymptotic of the solution. One should mention hcrcthat the term proportional to 1/r in (|12p gives a 
diverging contribution to the mass and is thus unphysical. Also, the constants ao, v?o must be chosen as 
zero, since we are dealing with perturbations. For solutions with <j = ao = <^o = 0, the asymptotic form 
(fl"2)) simplifies considerably and gives to the leading order : 
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Let us notice that the equations of perturbations, at each order in e and each massive modes are basically 
equations for waves propagating in the background of a, 6, / with a source terms. Let us write the 
perturbative expansion as X{r, z) = ^ i . Xf (r) cos(ikz)e J , X denoting the functions A, B, C. Then, for 

a given mode X\ , the equation is of the form 

nxf = j x} (x^<% (u) 

where □ is the d'Alembertian for the background and J x j is a polynomial of the X^. The exponent of each 
term in the polynomial times the order in e of the term must be equal to j in order to enter the equation at 
the e J order. 

The first mode containing a source term is the zero mode and decays as 1/r 5 . The mode X\ also decays 
as 1/r 5 , thus source terms for the order 3 in e will be of the form (Xl) 3 , X^X\ and decay as 1/r 15 . This 
decay is too quick to influence the leading order in the asymptotic expansion, thus leaving the free wave 
decay 1/r 5 . This argument can be extended by recurrence. 

4 Boundary Conditions and Numerical Techniques 

The background and first order boundary conditions are reminded here: 

a{r h ) = a , f(r h ) = , b{r h ) = , b'{r h ) = h, (15) 

where ao, b\ are real constants chosen such that the asymptotic solution fits the Feffer-Graham expansion 
©. For the first order fields, we have 

Ai(oo) = 0, <7x(oo) = 0, A' 1 (r h ) = A oo (A 10 ,C 10 ,a ,b 1 ) , 0[(r h ) = C^Aw, C 10 , a , b x ), (16) 

where Aoo, Coc arc the nonlinear conditions ensuring that the right hand side of ^ are regular at the 
horizon. The four conditions (|16|) arc not sufficient to fix the normalisation of A\ 1 C\. In order to solve this 
problem we add an extra equation expressing the constancy of the eigen value, (k 2 )' = 0, in order to solve 
the problem |13j and take advantage of the supplementary freedom to fix the arbitrary normalisation of F\ 
by means of Ci(rh) = 1. 

For the backreacting fields, the equations imply only one regularity condition for the equation of A' (rh), 
and the perturbations Aq, Co must vanish at infinity. 

A (oo) = 0, Co(oo) = 0, A^(r h ) = A (Co{rh),Aio,C 10 ,a ,bi), (17) 

where Aq is the regularity condition for Aq. This leaves a free parameter, for example C (rh), which is fixed 
imposing the physical decay (fH{|) . 

In practice, we integrate the background and first order with the solver Colsys, based on a Newton 
Raphson algorithm, with a mesh adapted in order to minimise the error |23j . The backreacting fields 
tends naturally to constants for large values of r. With the invariance ©, we give arbitrary values to 
Ao(r/j), Co(r/j), set A' (rh) according to the regularity condition and shoot using C' (rh) as a shooting 
parameter in order to fulfil the decay condition (|13p , transposed to the derivative of Co (r) . 

We used a Runge Kutta algorithm at order 4 to integrate the backreacting fields. The consistency of the 
result has been checked by comparing the result of the numerical integration with the analytic result of the 
asymptotic expansion. 

5 Properties of the solution 

5.1 Scaling properties and deformation parameter 

The equations present the following symmetry: 

I — * at , T}i — ► olT}i , k^k/a, (18) 
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for real values of a. Thus, given solutions with fixed AdS radius and varying horizon radius, one can convert 
these solution to fixed radius and varying AdS radius. This is usually referred to as copies of solutions. 
The scaling (|18[) affects the thcrmodynamical quantities in the following way : 

S^a 3 S , Th — > Th/ol , T — ► o?T , M — > a 2 M (19) 

Note that, in opposition to the asymptotically flat case, there is an intrinsic length scale in the model, 
given by I. As already mentioned in |13j . the solutions are characterised by two dimensionlcss quantities, 
/ii, fx 2 , with 

ML 
W = 167TGL--3 ' ^ 2 = 7 (20) 

with M the mass of the black string. 

Following |14j . one defines a deformation parameter, which quantifies how much the solution deviates 
from the uniform black string case. This parameter is given by: 

A= fimax_i (21) 
^min 

where i?max, -R m in represents the minimal and maximal areal radius at the horizon. Since we are dealing 
with perturbative non uniform black strings, our results will be valid for A << 1. The areal radius at the 
horizon is given by: 

R(z) = VmU r =w°M. (22) 



I T — Ty 

In the case of our interest, i?max = R{0) , -Rmax = -R(tt) which leads to: 

\ = 2C 10 e(l + eC 10 ) + O(e) 3 (23) 

5.2 Global charges and thermodynamical quantities 

The conserved charges of the solutions are computed using the counterterms approach [22]. This procedure, 
applied to uniform warped AdS black strings is explained in details in [12] , In this approach, one has to 
add suitable boundary terms S c t{h^ v ) to the action, in order to cancel the diverging terms, coming from the 
contribution AdS background. These terms do not contribute to the equations of motion and can be found 
in [22] . 

The variation of the total action Stot = S + S c t with respect to the boundary metric h^ v then leads to a 
regularized boundary stress tensor T^ u = — = . Then, a conserved charge 

Q ( = Jj a T ab n b dSi (24) 

can be associated with a closed four-surface S with surface element dS± and normal n b , provided £ a is a 
Killing vector of the boundary geometry. 

The first relevant corrections on thermodynamical quantities arises from the back reaction, since the 
integral (|24p involves an integration over z. However, terms with cos(fcz), cos(2fcz) factors cancel once 
integrated on z € [0, L], L = A straightforward computation gives the mass of the solution, which is the 
charge associated to the Killing vector d t : 

M = ^(-4c t + c z +30e\). (25) 

However, the solution's tension, which is the charge associated with d z is not affected by the perturbation, 
leaving the same value for the pNUBS than for the background: 

T=fgMc 2 + Ct) . (26) 
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Let us define also the relative tension n = TL/M. In the perturbativc approach, these quantities can be 
expanded in powers of e; i.e. keeping the first correction n + e 2 Sn , pi + e 2 (5/i;, i = 1,2, with the bar referring 
to the quantities evaluated on the background solution. After some algebra, we find the following relations 
involving the dimcnsionlcss quantities (|20[) : 

^ 1 = 30 £ 2 x , <5M2 = , ^ = _^I. (27) 

n [i x 

The entropy of the pNUBS, which is one quarter of the event horizon area, is given by 

c „ ^ T 2 3 (, (5ao (A 2 1Q - 12F 00 - 6A 1Q F 1Q - 9F 10 2 ) + 3 {A w - F 10 ) F 10 k 2 £ 2 ) A 

S = lyfolm r h ^1 — e j . (28) 

Depending on the sign of the correction, the pNUBS and thus presumably the full NUBS phase may be 
thcrmodynamically preferred rather than the uniform black string. 

The Hawking temperature as found by demanding regularity on the Euclidean section near the horizon, 
is given by 



r fl = ^(2* + 5^) (l + e ^ + ^ ^ - F " }F "' kV ) . (29) 

In 6 dimensions, the stress tensor has vanishing trace at the background level [12] , Let us finally mention 
that, as expected, this property holds for perturbativc non uniform black strings. 



6 Numerical results 

In the figures we present, we use the rescaling (fT8)) to choose the AdS radius to be £ — vTO without loosing 
generality. This value of £ corresponds to a unit negative cosmological constant. 

As already mentioned, the solutions are characterized by two dimensionless quantities, /ii, fj.2- In fig. [TJ 
we give the entropy per unit length as a function of the horizon temperature for the uniform black strings 
and the first relevant corrections for perturbative non uniform black strings, i.e. when the deformation 
parameter A << 1, for different values of /j,2- Since large black strings are known to be stable, we consider 
perturbative non uniform small black strings. Actually, the entropy per unit length of the background (the 
uniform phase) does not depend on the length of the black string. So this diagram is to be understood in 
the following way: each intersection of the non uniform black string phases with the uniform string case is 
the emanating point of a phase transition at fixed rescaled extradimension length [ii ■ 

For every cases we considered, the entropy of the non uniform phase is higher than the entropy of the 
uniform phase for the corresponding /X2 and the temperature decreases. For larger /X2, the entropy still 
increasing in the non uniform phase, but less than when [ii is small. This can be understood in the following 
way: since the length becomes of the order of £, so this gives "small - long black strings"; small for rh, large 
for L and the effect of the AdS radius is felt by the string in the z-direction. This is a consequence of the 
fact that the solution is characterised by two dimcnsionlcss constants. 
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Figure 1: The entropy per unit length as a function of the temperature for various (j,2- The entropy grows 
while the temperature decreases; for larger values of /i2, the entropy is nearly constant. 

Figure [2] is a plot of the relative tension n as a function of the rescaled mass fx\. This diagram is drawn 
in the same spirit as figure [TJ The relative tension of the non uniform phase is smaller than the uniform 
phase, but the mass is higher. 

As it should have been expected, the flat space picture seems to hold for AdS, at least for small-short 
pNUBS: non uniform black strings is preferred when m L. 




Figure 2: The relative tension as a function of the rescaled mass for various rescaled length. 
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7 Conclusion 



We constructed the first corrections to the thermodynamical quantities of AdS non uniform black strings in 
six dimensions. The main features of the phase diagram is similar to the asymptotically flat space case, except 
that the relative tension now depends on the rescaled mass for fixed rescaled length (n = n(/j,i), for fixed /K2). 

For the range of parameters we considered, the correction to the entropy was always positive, althought 
nearly vanishing for large rescaled length. The correction on the temperation was always negative, implying 
that the non uniform black string phase with small deformation parameter near the emanating point should 
be colder and more entropic, thus thermodynamically preferred. The tension-mass ration is decreasing in 
the non uniform phase, while the rescaled mass is increasing. It is remarkable that the tension is unaffected 
by the non uniformity, at least at second order in perturbation theory. 

It seems likely that small-short stable uniform black string should decay to caged AdS black holes, as it is 
the case in asymptotically flat space |17j while the unstable branch should decay to non uniform black string. 
We saw that the amount of change in entropy per unit length decreases when the rescaled length increases, 
but is still higher than the uniform phase. It was already known that the ratio r^jl plays an important role 
in the problem, but it seems that so does the rescaled length /z 2 . It seems likely that the string is affected 
by the AdS when the rescaled length is large i.e. when the length is of the same order than the AdS radius. 

Some interesting outlooks for this work should be to construct the full nonlinear solution for the warped 
AdS non uniform black string and the connection with the thermodynamics of the boundary CFT. Moreover, 
it would be interesting to generalise this work to d dimensions; in asymptotically flat space, there exists a 
critical dimension above which where the entropy of the non uniform branch increases instead of decreasing 
[8]. The consequence is that in cf < 13 the uniform branch is more stable than the non uniform branch. 
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